
METHOD FOR BISPECTRAL PICKING OF ANELLIPTICAL NMQ 
CORRECTION PARAMETERS 

The field of this invention is seismic prospecting. 
The invention relates particularly to processing of 
seismic traces in a common midpoint gather. 

The invention is more precisely related to a method 
of determining velocity V and anellipticity j\ parameters 
necessary to perform processing including NMO correction 
of seismic traces. 

Seismic prospecting usually consists of emitting 
seismic waves in the subsoil using one or several seismic 
sources, and making surface records ; of seismic data 
corresponding to seismic waves reflected on geological 
interfaces in the subsoil as a function of time (also 
called reflectors) using. receivers (also called geophones 
or hydrophones) depending on whether the prospecting is 
being done on land or at sea) and then processing these 
data to extract useful information about the geology of 
the subsoil. 

The record of the seismic energy produced by each 
receiver during data acquisition is called a seismic 
20 trace. 
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One conventional technique for seismic prospecting 
is multiple coverage, in which sources and receivers are 
arranged such that several seismic traces are grouped 
together at the same midpoint (in other words the point 
5 at equal distance between the source and the receiver for 
a given trace) . 

Although seismic traces contain useful information 
about seismic reflections and geology of. the subsoil, 
they also contain noise components. 
10 One of the first, purposes of the processing of 

seismic data is to eliminate, or at least attenuate, 
these unwanted noise components such that the useful 
information can be clearly identified and interpreted. 

One conventional method used to attenuate these 
15 noise components is Common Midpoint (CMP) gather. Traces 
with the same midpoint are then grouped as a function of 
the distance between the source and the receiver (called 
the offset) . 

In general, representation of seismic data in images 
20 requires the use of processing including: 

- a .so-called. TZO (Transform to Zero Offset) 
operation designed to compensate for the NMO effect of 
paths by correcting the arrival- time of reflections to 
the arrival time of traces with zero offset, 
25 - and a migration operation designed to reproduce, 

the correct shapes of the geological interfaces. 

Although these TZO and migration operations are 
usually carried out one after the- other, they can be 
carried out jointly. This is the case particularly when a 
3 0 time migration is done before stack (PSTM - Pre-Stack 
Time Migration) . 
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m a simpUfied menner, the Tz0 

r 1 ;" 00 °< Selsmic * —as and receivers 

placed at the common mid-point . 

5 oo'nt The P T Se " ^ aM reC ° rdS "^^-"ng the same 

rat " „ SUbSOil t0 inC " aSe the noise 
-- and the primary re(lectlQns ^ 

o I :° n ; ratl °' thUS benSfit tte Vantages 

or "multiple coverage". 

obUgueness correction method or « Worm*! M ove Out, 
correction method is use d. 

strj^ " ^ SUbSOU is Mentally 

■ s a , isd Wlth no laterai variation o£ y 

vlocxtx... lt can be shown that records that have the 

property of illumining the same point ^ ^ 

are records that have the same midpoint 

However, the image of a reflection in. the subsoil 
a™ at variable . times depending on the offset 
Therefore, in order to stack refactions, it is 'necessary 

o st art by correoting various recQrds ..^ ^ 

ail to a common zero offset reference, 

Historicslly, the NNO correction is based on a 
particuiarly simple model of the subsoil; a homogeneous 
model Wlth horizontal reflectors. 

25 In thiS m0de1 ' ^flections associated with each 

subsoil reflector are theoretical^ aligned along 
hyperbo as al so called indicetors, centred verticaliy 
above the midpoint. Y 

The arrival time of = ra r, . 
3 0 k •, • reflection is then a 

• JU hyperbolic function ,~,f *-v, 

runction of the • source-receiver offset, the 

shortest time being obtained at zero offset. ~ ~~~ 
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In order to . stack records in each gather, the NMO 

c» 

correction straightens the hyperbolas until they are 
theoretically horizontal. 

The NMO correction is then made based on the 
5 following hyperbolic equation for the travel time t after 
reflection, associated, with a source-receiver offset x 
pair-: 




in which t 0 represents the zero offset travel time and V 
10 denotes the average propagation velocity of waves in the 
subsoil . 

The simplistic model mentioned above is based 
particularly on assumptions such as low angles of 
incidence and an isotropic medium. 

15 But assumptions made for the simplistic model are 

too restrictive to describe a complex medium, and in 
particular cannot be applied to' the propagation of 
seismic Waves in an anisotropic medium (medium in which 
the - velocity of waves can vary, depending on the 

20 propagation direction)/' 

The use of a less simplistic model is necessary/ 
particularly due to: 

- the use of long streamers for acquisition in deep 
offshore., which iri particular means that traces with far 

25 offsets are acquired; 

- observation of anisotropy in clayey type 
sediments . 

The generally accepted assumption consists of 
modelling an anisotropic medium as a stack of 
3 0 transversely isotropic layers with a vertical axis of 



symmetry. This i s known as Vertical -rr a „ 

cal Tl ansverse Isotropy 
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(VTI) . 

Thus, it has been proposed to ^ 

f ^sea Co determine NMO 

corrections to be made by: 

5 - introducing vertical u 

a vertical non-homogeneity into a 

homogenous VTI medium TOd el, as presented ^ ^ 

I*™ 1 '" " ^ TS ~ »«. Velocity analysis 

for transversely isotropic media: G eo physi ce, 60 , 1550 . 

1 o 6 6 " ; 

10 or by 

- introducing VTI anisotrop y into . ^ ^ 

■ IL o 'I iS ° trOPiC 1SyerS ' " d — ^ document 
SIL1QI R. and BOUSQUIE H. , 20 oo, Anelliptic time 

processing based on an offset h^v-K i 

orrset hyperbola approach, 70 

Ann. Internet. „tg. : Soo . of Exp! . oeophys., 2245 . 224s „ 

TMs second approach combiaing vert . cal 
homogeneity and VTI anisotropy to give a new model of the 
subsoil, appears to be the best approach in mosc real 
.cases studied. 

The following equation for correction of the travel 
time t after reflection, ' using an anelliptic ■ ^ 
hyperbola, associated with a source receiver offset x 
pair, as derived from this model: 

t ( v,,„ = Jn_. + lfJIV~~? 

l + 8n VU+8nJ' (i+8,,)V 2 " Agnation (la) 
where V is the velocity conventionally used in seismic 

corresponding to small offsets snH „ ■ 

s ana 1 ls a parameter, 
called the anellipticity parameter. 

The document "SUAODEAU E. and SILIQI R . , 2001 
Anelliptic pre-stacK time migration. Annual International 
Meeting, CSEG Expanded Abstracts" a l S o proposes to 
include the NMO correction by anelliptic offset hyperbola 
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in the equation for the path time used during the PSTM 
migration operation. 

The PSTM migration equation is conventionally 
expressed in the form of a double anelliptic offset 
hyperbola, which is the sum of two square roots (the 
Double Square Root (DSQR) equation) . 

Taking account of anellipticity, the expression for 
this constant offset equation, becomes: 



t = 



8T) lft 0 / 2 Y (X-X m ^ /f t 0 / 2 V X-X m -h> 

! + 8ri 0 vl 1 + 8r lJ (1 + 8ti)V 2 vU + 8tiJ (U8r|)V 2 



10 Equation (lb) 

where : 

- the V and r\ parameters are those mentioned above, 
. - x m is the coordinates of the midpoints, 

- x - x m represents the migration aperture, 
15 - h is the half source -receiver offset, 

- t 0 is the double time for zero aperture of the 
operator. 

Note that when the aperture x - x m of the migration 

is zero, the PSTM correction equation (lb) becomes the 
20 NMO correction equation, (la). Thus, the NMO correction 

forms a special case of the PSTM migration; the case of 

PSTM migration with zero aperture. 

Therefore,, finally, in order to perform seismic data 

processing including NMO correction taking account of 
25 vertical heterogeneity and VTI type anisotropy, it is 

necessary to determine the velocity V and anellipticity r| 

parameters . 

' The said parameters V and T) can conventionally be 
estimated by carrying out two passes such that: 
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- during the first pass, the distribution of 
velocities V along the time axis is estimated, only using 
data with near offsets; 

- during the second pass, the anellipticity n is 
estimated along the time axis, using: 

- the distribution of velocities determined 
during the first pass, and 

- all data (including data with far offsets) . 
The document entitled: 

"SILIQI R . 2001, Technological leap in time 
processing focuses the data throughout anisotropic media: 
First Break, .19, No. 11, 612-618",' 

has also demonstrated that the parameters V and t) can be 
estimated in a single pass, during which bispectral 
15 analyses are carried out in order to simultaneously pick 

the V and n parameters along the time axis, using all 
data. 

However, dense analyses of NMO correction parameters 
are preferably made when the correction equation no 
longer depends on the time' t 0 (we then talk about static 
NMO correction) .. 

A static correction effectively provides a means of 
offsetting all samples each forming traces of the same 
time 5t, for a given offset. 
25 Thus, when a static correction is made, the number 

of calculations to be made can be significantly reduced 
and the trace stretching phenomenon is eliminated, which 
makes the said dense analyses viable. 

At the moment, the only way to make velocity and 
30 anellipticity analyses is to use two passes, in 
particular making use of parabolic approximations of NMO 
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residuals, in order to obtain a dense picking of the V 
and r) parameters. 

In this context: 

- the first step is to estimate residual velocities, 
using near offset data, according to a first estimate of 
velocities; 

- the next step is to estimate anellipticity on all 
data, using velocity updates done previously. 

A mute law must also be defined to estimate residual 
velocities, so as- to only keep data that are considered 
as being near offset data. 

However, the efficiency of the two-pass analysis is 
particularly sensitive to the choice of such a mute law. 

On the other hand, the anellipticity tj is estimated 
mainly making use of far offset data. 

But, the anellipticity estimate made for far offset 
data is not very precise, such that the correction made 
is finally inaccurate. 

One purpose- of- the. invention . is . to overcome these 
limitations and disadvantages, by proposing a processing 
including a static NMO correction which is more efficient 
and more precise. 

More precisely, the purpose of the invention is a 
dense determination of velocity V and anellipticity r\ 
parameters in a single pass using all available data, in 
other words making use of all offset ranges. 

To achieve this, the invention proposes a method of 
determining the velocity V and anellipticity r\ parameters 
for processing seismic traces from a common midpoint 
(CMP) gather including an anelliptic NMO correction, 
comprising: 



- a preliminary step to define a plurality of nodes 
(dtn, x 0 ) , the said nodes being indicative of parameters 
dtn and x 0 representing the NMO correction for the 
maximum offset and the zero offset travel time in 
hyperbolic coordinates, the said preliminary step being 
followed by 

- for each node (dtn, x 0 ) defined in the preliminary 
step, steps for the NMO correction of traces in the CMP 
gather as a function of the values of the said parameters 
dtn and x 0 at the. node considered, and for calculating 
the semblance . function associated with the said NMO 
correction for the node considered; and 

- for each picked time t 0 , a step including 
determination of the maximum semblance node (dtn(t 0 ), 
to (t 0 )), 

- and a final step to convert the dtn(t 0 ) and t 0 (t 0 ) 
parameters so as' to obtain the velocity V(t 0 ) and 
anellepticity ii (t 0 ) laws. \ 

According to a. first- embodiment of the invention, 
the processing 'done is a' static NMO correction of seismic 
traces. 

According to a second embodiment of the invention, 
the processing done is a PSTM migration of seismic 
traces, the said PSTM migration including a static PSTM 
NMO correction of the said seismic traces. 

A preferred but non- limitative aspect of the process 
according to the invention relates to the definition of 
the dtn and x 0 parameters relative to the velocity V and 
the anellipticity r| so as to make static NMO corrections, 
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according to dtn = Sti t + 
to 

1 + 8ti ' 





2 X > 




(3+8ri)V 2 



and t 0 = 



Other aspects, purposes and advantages of this 
invention will become clearer after reading the following 
detailed description with reference to the attached 
Figures in which: 

- Figure la shows the anelliptic offset hyperbola 
used to make the NMO correction and illustrates s the 
meaning of the parameters x 0 and dtn; 

- Figure lb represents the DSQR equation of the 
double anelliptic offset hyperbola of the PSTM migration 
and illustrates the meaning of the t 0 and dtn parameters ; 

- Figure 2 illustrates the effect of the t 0 
parameter of the reflection curve corrected by the NMO 

15 correction; 

- Figures 3a and 3b represent the analysis volume 
(t 0( dtn, T 0 ) in which the bispectral picking of the dtn 

• and t 0 parameters according to the invention' is made; 

- Figure 4 combines each of (V, Van ) and (dtn, x 0 ) 
approaches, representing their corresponding bispectral 
analysis panel; 

- Figure 5 shows the correspondence, between the 
(dtn, t 0 ) and (V, V^) pairs of NMO correction parameters; 

- Figure 6 shows a CMP gather of true seismic traces 
before the NMO correction and the bispectral picking of 
the (dtn, T 0 ) parameters corresponding to this gather of 
traces; 
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- Figure 7 shows the semblance function and velocity 
and anellipticity functions deduced from the picking of 
the dtn and x 0 parameters in Figure 6. 

- Figure 8a is a flowchart representing the steps of 
a first particular embodiment of the invention, namely 
determination of the V and r| parameters to make an NMO 
correction; 

- Figure 8b is a flowchart representing the steps of 
a second particular embodiment of the invention, namely 
determination of V- and r) parameters to do a PSTM 
migration; 

- Figure 9 is a diagram illustrating the different 
operations carried out in order to determine the V and r\ 
parameters to make a PSTM migration. 

In general, the method according to the invention is 
a method for processing records of variable offset 
seismic traces, and this processing uses the recorded 
seismic traces to build common midpoint (CMP) trace 
gathers, and the traces in each gather are subjected to 
an NMO. correction. 

In particular, the method according to the invention 
determines the velocity V and anellipticity r| parameters 
to perform processing including such an NMO correction of 
the seismic traces of a CMP gather. 

For example, the said processing may be: 

- an anelliptic NMO correction of the seismic 
traces'; 

- an anelliptic PSTM migration which, as we have 
seen above, jointly applies TZO and migration operations 
(in the following we will refer to a PSTM NMO 
correction) . 
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The following description is more specifically 
applicable to the NMO correction. However, considering in 
particular the two particular embodiments of the 
invention that will be described below, we can see that 
this description is equally applicable to any processing 
including an NMO correction, particularly processing 
including a PSTM NMO correction. 

Two new parameters dtn and x 0 are considered in 
order to specifically determine V and rj : 

- T 0 that represents the zero offset travel time in 
"hyperbolic coordinates" (see Figure 1) , 

to = ■ ~ T| Equation (2) 

l + 8r) 

and, 

- dtn that represents the NMO correction for the 
largest offset (see Figure 1) , 

dtn = t x=x „„ - t x « 0 namely: 




x. 2 



dtn = 77c ° + J 7~T" + "^-T Equation (3) 

It is important to note that dtn is defined relative 
to the velocity V and the anellipticity r\, while x 0 is a 
perfectly anelliptical parameter defined relative to the 
anellipticity rj, independent of V. 

Accordingly, the velocity V and anellipticity n. 
parameters may be calculated in accordance with equations 
(2) and (3), using the following conversion equations (4) 
and (5) : 



y __ max 



Equation (4) 



idtn(dtn + 2T 0 )^ 
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and ii = - Equation (5) 



Figure la shows the anelliptical offset hyperbola 
used to make the NMO correction and illustrates the 
meaning of the x 0 and dtn parameters. 

The said "hyperbolic coordinates" are shown in this 
Figure 1. Their origin on the time axis is taken at the 
intersection of the said time axis with the asymptote 
tangent to the said far offset offset hyperbola (see 
equation (la) ) . 

Using the (dtn, x 0 ) parameters,, equation (la) for 
the anelliptical offset hyperbola becomes: 

. „ (I dtn(dtn + 2t n ) , 

t - to - To + K+ i r-^x 2 Equation (6a) 

V x max 

15 The ( dtn > ^o) parameters defined for the velocity V 

and. anellipticity v\ are thus used to make the NMO 
correction CORIW = t - t 0 to be : applied to offset x 
traces independent of t 0 : 
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COR Rnmo (x) =—c 0 + T 2 + dtn(dtn + 2T 0 ) x2 Equation (?a) 



X max 



Therefore, this is a static NMO correction. In other 
words, the data recorded on a given offset trace will all 
be corrected in the. same way for a (dtn, r 0 ) pair, 
independently of the time at which these data were 
acquired. 

2 5 Consequently, the estimate of the velocity and 

anellipticity parameters is not disturbed by stretching 
of traces usually observed when dynamic NMO corrections 
are made. 
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Similarly using the two parameters (dtn, t 0 ) , 
equation (lb) for the DSQR double offset hyperbola for 
the anelliptic PSTM migration becomes: 



5 t - t. - t. + N , dt "(dtn + 2T n )(x- Xm + by 

V 4 x? """ * 



v max 



/T^ + dtn(dtn + 2T 0 )(x-x m -h> 
V4 x 2 

max 



|2 

Equation (6b) 



Figure lb shows the DSQR double offset hyperbola of 
the PSTM migration defined with the (dtn, x 0 ) parameters. 

In the context of the PSTM migration, x majc represents 
the maximum offset and aperture of' the migration. 

The (dtn, t 0 ) parameters defined relative to the 
velocity V and anellipticity r) are thus used to make the 
PSTM NMO correction CORR PSTM = t - t 0 to be applied to the 
15 offset x traces independent of t 0 (equation 7b) ) : 

Corr PSTM (x) = -T 0 + . R " dtn(dtn + 2Tt){x-x+ hy + \*l , dtn{dtn + 1t 0 ){x - x - hf 

" ' V 4 ~lL ~ V4 7 

Therefore,' this is a- static PSTM correction. In 

other words, the se of aperture samples x - x m of an 

"iso-offsef cube are offset by the same time for a given 
20 (dtn, T 0 ) pair. 

By analysing the (dtn, t 0 ) parameters in' several 
picking times, the method' according to the invention can 
be used in particular to determine the (V, r,) parameters 
necessary for processing including an anelliptic NMO 
25 correction of the traces of a CMP gather. 

The said process includes the steps presented below 
in a simplified manner. 
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During the preliminary step, an analysis volume is 

defined including several nodes (dtn, t 0 ) . 

The following are performed for all nodes in this 
volume : 

- firstly, according to equation (7a), the static 
NMO correction of . traces in the studied CMP gather, as a 
function of the values of the dtn, t 0 parameters at the 
node considered, the said static correction being valid 
for any picking time; 

- secondly, the semblance as a function of time 
associated with the correction made in the previous step 
is calculated. 

Finally, the following are determined for each time 
in the several picking times: 

- firstly, the node (dtn, x 0 ) used to make an 
optimum correction, for example with regard to the 
semblance criterion (this type of criterion typically 
being used in seismic processing" ' to • "measure the 
horizontality" of reflection curves and to determine the 
reliability of the. picking) ; ■ - 

- secondly, the values of the dtn, t„ parameters at 
the said maximum semblance node are converted to values 
of the velocity V and anellipticity r\ parameters at the 
said picking time considered. 

The velocity law (in other words all (picking time, 
V)) pairs and the anellipticity law (all (picking time, 
r))) pairs are thus - defined. 

Finally, the NMO correction for all seismic traces 
may be made using, these velocity V and anellipticity r] 
laws in equation- (la) for the anelliptic offset 
hyperbola. 
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Similarly, and as will be described in detail later, 
the (dtn, x 0 ) parameters of the static PSTN migration 
(see equation 7b) ) may also be checked. The velocity V 
and anellipticity 7] laws are then determined and may be 
5 used to make the PSTM migration in the DSQR equation (lb) 
for the double anelliptic offset hyperbola. 

Figure 2 illustrates the effect of the parameter t 0 
(and therefore the anellipticity r| according to equation 
(2) on the curvature residues after the NMO correction. 
10 Note that the vertical time scale of the curve in 

Figure 2 is exaggerated so that this effect can be 
clearly understood. 

Three curves are shown in Figure 2, for which the 
dtn parameter is fixed to the correct value and the 
15 parameter x 0 is assigned different values. 

The central curve shows the case in which t 0 is 
equal, to its. correct value t 0 i, in other words when the 
corresponding, anellipticity t)i is equal to its true-value 
"Htrue- As expected, the corrected' reflection curve is then 
2 0 horizontal.. 

The uppexv curver shows, the case in which t 0 : is equal 
to a value x Q2 . less than its correct value x 0 i, the 
corresponding anellipticity- T) 2 being greater- than its 
true value r) t rue- 
2 5 The lower curve represents the case in which t 0 is 

equal to a value x 03 greater than its true value x 0 i, the 
anellipticity r|3 being less than its true value T^me- 

Note that these lower and upper curves show that the 
"horizontality" of the corrected reflection curve is 



acceptable at near offset (x « 0) and at far offset (x * 

On the other hand, significant curvature residues 
are observed when the offset x no longer tends towards 
one of these limiting values 0 and x^. In particular, 
particularly significant residues are observed for an 
offset x centred in the middle of the offsets range. 

For example, and as is shown in Figure 2, when t 0 is 
equal to t 02 , an RMOt 02 residual correction must be made 
to the offset traces x RMO To2. Similarly, when t 0 is equal 
to t 03 , a residual correction RMOx 03 must be made to 
offset traces XrmoT 03 . 

Due to these significant curvature residues, the 
offsets range can be used almost in its entirety in order 
to determine the anellipticity tj. 

Therefore, setting parameters in (dtn, x 0 ) for the 
NMD correction make it possible to use available data, for 
all offsets (x between 0 and x max ) in determining the 
anellipticity T|. '"' ' "• ■ "* 

As has already been mentidnedy this" is not": the "case 
for NMO corrections for which parameters have.;' been-- set' 
with V and r| for which the estimate of the anellipticity' 
T) is made essentially using far offset data. 

Thus, the effect of the new anelliptic parameter t 0 
is distributed on all offsets, unlike the anellipticity r| 
that only affects far offsets. Therefore, the "behaviour" 
of the parameter t 0 gives a better constraint on 
anellipticity values. 
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As already mentioned above, the optimum dtn and t 0 
parameters are determined within a 3D analysis volume 
(t Q , dtn, t 0 ) . 

Several nodes (dtn, t q ) , in other words several 
pairs of dtn, x 0 parameters for which values are known, 
are considered in the said analysis volume. 

Nodes are usually regularly spaced from each other, 
by an increment Adtn on the dtn axis and an increment At 0 
on the t 0 axis. 

Minimum values dtn min , T 0min , t 0mi „ and maximum values 
dtn max , T 0ma x/ t 0 max of the To and to respectively provide a 
means of defining the limits of the said analysis volume. 

Advantageously, plausible values of the velocity 
parameter V and the anellipticity parameter tj may be used 
to define a corridor [dtn mia (t 0 ) , dt nmax (t 0 )], [T 0min (t 0 ), 
tomax(to)] inside the said analysis volume. 

This corridor restricts the analysis volume and 
therefore the number of nodes (dtn, t 0 ) that have to be 
considered to determine the optimum pair (dtn, t 0 ) '. 

Although use of the said corridor . is beneficial for 
the efficiency of the process according to the invention, 
it also provides a means of constraining solution, towards 
the right phenomena, without needing to consider 
incompatible (dtn, t 0 ) (and therefore V, i\) ) pairs, for 
example related to multiple reflections or miscellaneous 
interference phenomena. 

Figures 3a and 3b represent the said analysis volume 
(t 0 , dtn, T 0 ) in the context of an example of true seismic 
data processing during an NMO correction done according 
to the process according to the invention. 



19 



Figure 3a represents three 2D panels a, b and c of 
the analysis volume: 

- panel a is a panel (dtn, t 0 ) with constant t„; 

- panel b is a panel (dtn, t 0 ) with constant x 0 ; 

- panel c is a panel (x 0> t 0 ) with constant dtn. 
Figure 3b diagrammatic ally shows the 3D analysis 

volume (t 0 , dtn, x„) and three intersections with this 
volume along three planes with constant t 0 , x 0 , dtn 
respectively, each of these intersections being projected 
onto the corresponding side a, b or c in Figure 3a. 

Panels (dtn, x 0 ) with constant t 0 (panel a in the 
above example) are panels in which the bispectral picking 
of the dtn and x 0 parameters is made, for example 
according to the maximum NMO correction semblance 
criterion, for the considered picking time t 0 . 

The [dtiWt 0 ), dtWt.)] [tomin (t 0 ) W (t 0 J ] 
corridor mentioned above for the effective analysis 
inside the analysis volume is also shown in Figure 3a 
(lighter sone on each panel) . 

The seismic data sampling time step At 0 defines the 
difference between two successive picking times (times 
for which the maximum associated semblance node (dtn, x 0 ) 
is determined) and therefore the number of bispectral 
picking panels to be considered. 

Automatic picking also provides a means of 
extracting the dtn, x 0 parameters at a density that is 

greater when the At 0 increment between picking times is 
small . 

Sampling of the analysis parameters dtn and x 0 is 
directly related to the resolution of the seismic ' 
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exploration; dtn and x 0 effectively have the same 
dimensions as seismic records. 

The systematic search for the maximum semblance, 
conventionally known in itself, provides a means of 
5 determining the (dtn, t 0 ) pair providing best focus, for 
a given picking time t 0 . 

Parabolic interpolations about the values of nodes 
(dtn, t 0 ) can also provide a means of evaluating values 
of the dtn, t 0 parameters between the different nodes 
10 that have actually been picked. And in particular this 
type of evaluation makes determination even more precise 
(in contrast to determination limited to nodes in the 
corridor) of the dtn, t 0 parameter pair maximising the 
semblance function. 
15 Finally, the velocity V and anellipticity rj 

parameters are determined, always for the picking time t 0 
considered, by using, the above mentioned conversion 
equations (4) and. (5). • 

Figure 4 shows the comparison between the two 
20 approaches" (V, VaJ and (dtn,,.T 0 ), by representing the 
bispectral picking panel for each, for a given picking 
time. 

The figure at the right illustrates the conventional 
approach (V, V an ) for which the two axes are the velocity 
25 axes (the anellipticity r| being related to the ratio of 



ifv 4 N 

these two velocities according to n = — -1 

*U 4 . 



) 



The left figure illustrates the approach (dtn, t 0 ) 
according to the invention for which the two axes are 
time axes. 



It is important to note from the study in Figure 4 
that the dtn and x Q parameters appear to be decorrelated. 
This "decorrelation" is striking when the two approaches 
are compared, spreading of the (dtn, w x 0 ) spectrum 
actually being much narrower than spreading of the (V, 
Van) spectrum. 

Therefore the picking made in the context of the 
(dtn, t 0 ) approach according to the invention is more 
precise than the approach carried out in conventional 
methods . 

Furthermore, this decorrelation enables filtering of 
the dtn and x 0 pickings separately, while keeping NMO 
corrections. This is not the case for the V and r\ 
parameters for which a reduction in one of the parameters 
must necessarily be compensated by an increase in the 
other and vice versa. 

And due to interpolations and individual filtering 
of parameters according to the invention, dtn and t 0 , it 
is then possible to perform simultaneous, interpolation 
and filtering' of the standard NMO correction parameters V 
and T| . : 

Figure 5 shows' non-linear correspondence-"" according 
to- equations (4) and (5) mentioned- above', between the 
pair of time parameters (dtn,. t 0 ), and the pair. of. 
velocity parameters (V, V an ) . 

Figure 6 shows, from left to right: 

- a CMP gather of real seismic * traces before NMO 
correction ; 

- the picking of the dtn parameter corresponding to 
this gather of traces; 
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- the picking of the t 0 parameter corresponding to 
this gather of traces. 

The straight line x 0 = t 0 on the picking of the 
parameter t 0 at the right in Figure 6, corresponds to the 
purely hyperbolic reflection curves. 

From right to left, Figure 7 represents the 
semblance function and the velocity V and the 
anellipticity r\ functions (see equation (4) and (5) ) 
deduced from the' picking of the dtn and t 0 parameters 
shown in Figure 6 . 

On this example of real seismic data processing, it 
can be seen that the values of V and T] obtained generally 
corresponding to a semblance of more than 4 0%. 

The following description contains details of two 
particular embodiments of the invention. 

The first of these modes relates to a process for 
determining optimum parameters to make an anelliptic NMO 
correction to the traces of a CMP gather' (see the 
different steps- shown in the= flowchart in Figure 8a) . : 

With reference- to Figure 8a, this first embodiment 
• includes an initialisation- step la during which the 
following- operations are carried out in sequence: 

- determination of the limits of the analysis volume- 

[dtn m i n/ dtn maX ] [^Ornin, T 0 max] t^Omin/ t 0 max] } 

- calculation of NMO corrections CORR NMD (equation 
(7a)) for all offsets and for all nodes (dtn, x 0 ) 
included in the analysis volume; 

- delimitation inside the analysis volume of the 
corridor [dtn min (t 0 > , ' dtrw (t 0 > ] [T 0min (t 0 ) w (t 0 ) ] of 
plausible velocity and anellipticity values. 
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Once the initialisation step la has been carried 
out, a step 2a to calculate velocity V(t 0 ) and 
anellipticity x] (t 0 ) laws is performed for each gather of 
CMP traces. 
5 This step 2a comprises: 

- a first operation 3a carried out for each node 
(dtn, t 0 ) in the corridor defined in the initialisation 
step la, during which the following operations are 
carried out in sequence for each picking time- t 0 : 
10 ' application of static NMO corrections CORR,^ 

for all offsets along the corridor, pr'ecalculated during 
the initialisation step la; 

- calculation of the semblance function on data 
corrected along the corridor using a time window 

15 appropriate to the dominant wavelet; 

- summation ("stack" calculation) of data 
corrected along the corridor (only near offset data can 
be used advantageously for this purpose); • . 

- a second^ operation 4a carried out for each picking 
time t 0 (the said times being at intervals of At 0 between 
[tomiiw. tomaxJ during,, which-, the ■=' following, operations are 
carried out : ... 

- search for maximum semblance, in the [dtn min (t 0 ) / 
dtn max (t 0 )]. [T 0m in (t 0 ) T 0max (t 0 ) ]. corridor of . the bispectral 

25 panel (dtn, t 0 ) ; 

- check the fact that the position in (dtn, t 0 ) at 
the maximum semblance corresponds to a summation extreme 
value for the same values dtn and t 0 ; 

- creation of the dtn(t 0 ), t 0 (t 0 ) and semblance 
30 ' (t 0 ) series; 



- a third operation 5a designed to select and adjust 
the pickings obtained, during which the following 
operations are carried out: 

- increasing sort of the semblance series (t 0 ) ; 
validation of pickings dtn and x 0 for which the 

distance in time to the highest semblance pickings is 
greater than a predefined value; 

- adjustment of picking dtn and x 0 values 
validated by parabolic interpolations using surrounding 
values; 

- retention of picked, validated and adjusted 
values if it is possible to calculate the Dix interval 
velocities with the pickings with the highest 
similarities . 

- a fourth operation 6a designed to use equations 
(3) and (4) to convert picked, validated and adjusted and 
retained values of dtn and x 0 during operation- 5a, into 
velocity V and anellipticity r| laws. 

The velocity V" ' and anellipticity r\ laws as a 
function - of- time are thus- perfectly determined. The 
ahelliptical NMO r correction of seismic traces in the CMP 
gather can thus be done precisely. 

The second' particular embodiment of the invention 
relates to a process for determination of optimum, 
parameters for anelliptical PSTM migration of traces in a 
CMP gather. 

This second embodiment may be included in a 
generalisation of the first mode discussed above. 

As has already been demonstrated, the use of the 
(dtn, x 0 ) parameters, enables static PSTM NMO corrections 
(see equation (7b)). 
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In the context of a PSTM NMO correction, this 
application has the same advantages as previously 
discussed for the NMO obliqueness correction. 

More precisely, it will be noted that the first 
5 embodiment is ■ only a special case of the second 
embodiment corresponding to the case of a zero migration 
aperture . 

With reference to Fi gure 8b , the second embodiment 
-eludes an initialisation step lb during which the 
3 following operations are carried out in sequence: 

- determine the • limits of ' the analysis volume 
[dtn minf dtn max ] , [tomia( w] # [tomia# tomax] ; 

- calculate NMO corrections CORR PSTM (equation (7b)) 
for all nodes (dtn, x 0 ) inc l u ded in the analysis volume 
and for all migration offsets inside the migration 
aperture; 

- delimitation of plaU sible velocity an<J 
anellipttcity values inside the analysis voiume of the 
corridor [dtn^ , w ... dtn TO ,t 0 ) , f w , w . ^ (t§jj . , 
. Once this initialisation step ih has been done, the 
sard first embodiment creates a step 2b to calculate the 
velocity v«t,) and- anellipticity „ (ti) . law3 for ^ 
gather of CMP traces.. 

This step 2b comprises: 

- a first operation 3b done for each node (dtn, x 0 ) 
of the corridor defined during the initialisation step 
la,, during which the following are carried out in 
sequence : 

- for each offset class, operations for: 

- application on all midpoints inside the 
migration aperture along the corridor, of static 
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corrections C0RR PSTM precalculated during the 
initialisation step lb; 

- summation of corrected midpoints along the 

corridor. 

- for each picking time t 0 , operations to: 
. - calculate the semblance function on 

corrected data along the corridor using a time window 
appropriate to the dominant wavelet; 

- summation (stack ■ calculation) of data 
corrected along the. .. corridor (only data with small 
offsets can be used advantageously for this purpose) : 

- a second operation 4b carried out for each picking 
time t 0 , similar to. the operation 4a described above, to 
create the dtn(t 0 ), x 0 (t 0 ) and semblance t 0 series; 

15 ' a third operation similar to operation 5b 

described above, designed to select and adjust the 
pickings obtained; 

- a fourth operation; 6b similar' to operation 5a 
described above,, designed to convert- values, of dtn, t 0 

20 into velocity. V and anellipticity r| laws.. 

The velocity Vi, and anellipticity n laws as a 
function of . time are thus, perfectly determined. And the 
PSTM migration of seismic traces in the CMP gather may 
thus be done accurately. 

Figure 9 illustrates the second embodiment of the 
invention that has just been described. 

Seismic data are initially grouped into iso-offset 
cubes . 

As has just been described, NMO corrections CORR PSXM 
are applied (see first operation in step 3b) for each 
offset class (namely for each iso-offset cube) on all 
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midpoints. The arrow marked with the MIG label in Figure 
9 illustrates this operation. 

The midpoints thus corrected are then summated 
during a second operation in step 3b, the arrow labelled 
5 STACK in Figure 9 illustrating this operation. 

These two MIG and STACK operations are specific to 
the second embodiment of the invention (PSTM migration) . 

The following operations are performed in the 
context of each of the two embodiments discussed (the NMO 
10 obliqueness correction, as we. have described, 
corresponding to the special case of a PSTM migration 
with zero migration aperture) . 

The semblance is then calculated for each picking 
time (second operation in step 3a for the NMO correction 
only, third operation in step 3b for the PSTM migration) , 
the arrow labelled "semblance" in Figure 9 illustrating 
this operation. 

Obviously, the MIG, STACK and semblance calculation 
operations are implemented for each node (dtn,' x 0 ) . 

The "automatic bispectral picking" mentioned in. 
Figure 9 corresponds to picking of maximum semblance 
parameters (dtn, x 0 ) for each picking time t 0 (operations - 
4a, 4b in Figures 8a, 8b respectively)'. ' The arrow 
labelled MAX illustrates the search for the maximum 
25 semblance. 

Finally, the pickings of the picked parameters (dtn, 
to) are converted into velocity V and anellipticity r\ . 
laws (operations 6a, 6b in Figures 8a, 8b and arrow 
labelled CONV in Figure 9) . 
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